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Brief introduction to the discrete quantum mechanics is given together with the 
main results on various exactly solvable systems. Namely, the intertwining rela- 
tions, shape invariance, Heisenberg operator solutions, annihilation/creation op- 
erators, dynamical symmetry algebras including the g-oscillator algebra and the 
Askey- Wilson algebra. A simple recipe to construct exactly and quasi-exactly solv- 
able Hamiltonians in one-dimensional 'discrete' quantum mechanics is presented. It 
reproduces all the known ones whose eigenfunctions consist of the Askey scheme of 
hypergeometric orthogonal polynomials of a continuous or a discrete variable. Sev- 
eral new exactly and quasi-exactly solvable ones are constructed. The sinusoidal 
coordinate plays an essential role. 
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1. Introduction 

This is an expanded version of my talk in Robin BuUough Memorial Symposium, 11 
June 2009, Manchester 1 met Robin BuUough in the Niels Bohr Institute, Copen- 
hagen in 1979, after finishing two papers on unified geometric theory of soliton equa- 
tions (Sasaki 1979a, b). Robin was my mentor and navigator when I started cruising 
into an unknown sea of integrable systems and soliton theory. Thirty years ago, the 
subject was young and rough but it was very fertile and full of good promises. He 
gave me good suggestions and comments on my third and fourth papers on geomet- 
ric approach to soliton equations (Sasaki 1979c, 1980a). We produced joint papers, 
one on 'Geometry of AKNS-ZS scheme' (Sasaki & Bulloughl980) and the other 
on non-local conservation laws of sine-Gordon equation (Sasaki & BuUough 1981). 
Collaboration was a fruitful experience and I learned a lot from him and got to 
know many friends and experts through him. As a result, I am still working in the 
general area of integrable systems, and in particular, the general mechanism of in- 
tegrability and solvability as a means to understand nonlinear interactions beyond 
perturbation, rather than some specific properties of certain equations. 

f Contribution to the Theme issue of the Philosophical Transactions A entitled Nonlinear 
phenomena, Optical and Quantum solitons, dedicated to the figure of Robin BuUough. 
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Here I present an overview of the problems I have been working in recent twelve 
years; exactly and quasi-exactly solvable quantum particle dynamics. Due to the 
lack of space, I have to concentrate on the most basic but most interesting part of 
the problems; the quantum mechanics (QM) of one degree of freedom. The subjects 
to be covered in this article include: general setting of discrete quantum mechanics 
with pure imaginary shifts and real shifts, shape invariance, exact solvability in 
the Heisenberg picture, dynamical symmetry algebras, unified theory of exactly 
and quasi-exactly solvable discrete quantum mechanics. Let me emphasise that 
the basic structure of integrability and methods of solutions are shared by the 
ordinary quantum mechanics (oQM) and the discrete quantum mechanics (dQM). 
For each specific subject, I will present basic formulas common to the ordinary and 
the discrete QM and then provide explicit forms of various quantities which are 
different in the oQM or the two types of dQM. 

One underlying motivation of this series of works is to inaugurate a new the- 
ory for the Askey scheme of hypergeometric orthogonal polynomials (Andrews et 
al 1999; Nikiforov et al 1991; Ismail 2005; Koekoek & Swarttouw 1996) within 
the framework of quantum mechanics, so that its abundant concepts and methods 
would be available for the research of orthogonal polynomials. So far this project 
has been extremely successful. Every Hamiltonian system describing the hypergeo- 
metric orthogonal polynomial is shown to be shape invariant thus exactly solvable 
in the Schrddinger picture (Odake & Sasaki 2005a, 2008b, c). As a byproduct a 
universal Rodrigues type formula (j2.45l) for all the hypergeometric orthogonal poly- 
nomials is obtained. It also turned out that all these Hamiltonian systems describing 
the hypergeometric orthogonal polynomials are exactly solvable in the Heisenberg 
picture (Odake & Sasaki 2006a, b). Namely, the Heisenberg equation for the si- 
nusoidal coordinate ri{x) is exactly solvable. The positive/negative energy parts of 
the exact Heisenberg solution for rj{x) provide the annihilation/ creation operators^ 
which together with the Hamiltonian constitute the dynamical symmetry algebra 
of the exactly solvable system. The q-oscillator algebra is obtained as a dynamical 
symmetry algebra of the system for the q-Hermite polynomials (Odake & Sasaki 
2008a). The relationship dictated by the annihilation/creation operators can be 
interpreted as quantum mechanical disguise of the three term recurrence relations 
of the orthogonal polynomials. Various coherent states as the eigenvectors of the 
annihilation operators are explicitly obtained (Odake & Sasaki 2006a, b, 2008c). 
This project also opens a rich frontier for various deformed polynomials in terms 
of modification of Crum's theorem in the ordinary QM (Krein 1957; Adler 1994) 
and for the discrete QM (Odake & Sasaki 2010b). The problem of determining 
the deformed weight functions can be solved without recourse to the cumbersome 
moment problems. One recent highlight of the project is the discovery of infinitely 
many exceptional orthogonal polynomials related to the Laguerre, Jacobi, contin- 
uous Hahn, Wilson and Askey- Wilson polynomials (Odake & Sasaki 2009c, d, e, 
2010a, b; Ho, Odake & Sasaki 2009). 



2. 'Discrete' Quantum Mechanics 

Throughout this paper we consider one-dimensional quantum mechanics. The dy- 
namical variables are the coordinate x and its conjugate momentum p, which is 
realised as a differential operator p — —i-^ = —idx- The discrete quantum mechan- 
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ics is a generalisation of quantum mechanics in which the Schrodinger equation is 
a difference equation instead of differential in the ordinary QM (Odake & Sasaki 
2008b, c). In other words, the Hamiltonian contains the momentum operator in 
exponentiated forms e^^^ which work as shift operators on the wavefunction 

e'^'^Pilj{x) =ij{xTil3). (2.1) 

According to the two choices of the parameter P, either real or pure imaginary, 
we have two types of dQM; with (i) pure imaginary shifts /3 = 7 G M^^o (pdQM), 
or (ii) real shifts = i (rdQM), respectively. In the case of pdQM, ^p{x =F *7), we 
require the wavefunction to be an analytic function of x with its domain including 
the real axis or a part of it on which the dynamical variable x is defined. For the 
rdQM, the difference equation gives constraints on wavefunctions only on equally 
spaced lattice points. Then wc choose, after proper rcscaling, the variable x to be 
a non-negative integer, with the total number either finite {N + 1) or infinite. 

Here we will consider the Hamiltonians having a finite or semi-infinite number 
of discrete energy levels only: 

= 5(0) <£(!)< £(2) <••• . (2.2) 

The additive constant of the Hamiltonian is so chosen that the groundstatc energy 
vanishes. That is, the Hamiltonian is positive semi-definite. It is a well known 
theorem in linear algebra that any positive semi-definite hermitian matrix can be 
factorised as a product of a certain matrix, say A, and its hermitian conjugate ■ 
As we will see shortly, the Hamiltonians of one-dimensional quantum mechanics 
have the same property, either the ordinary QM or the discrete QM both with the 
imaginary and real shifts. 



(a) Hamiltonian 

The Hamiltonian of one dimensional quantum mechanics has a simple form 

H = A^A, (2.3) 

in which the operators A and A^ are: 

^ , def d dwix) .+ d dw(x) , ^ ^ ,^ 

oQM A=- A^ = -- w(x)gR, (2.4) 

dx dx dx dx 

pdQM A =*i(eTf/V^*(a;) - e-'P/'^^V{x)), V{x),V*{x) e C, 

A'' = -i{^/V(xje'^P/^ -^/W(xje-''P/^), (2.5) 

rdQM A = ^/B{x) - e^^D{x), B{x) > 0, D{x) > 0, 

A^ = ^B{x) - ^D{x) e-^, D{0) = 0, B{N)=0. (2.6) 

Here V*(x) is an analytic function of x obtained from V{x) by the *-operation, 
which is defined as follows. If f{x) = ^a„a;", a„ G C, then f*{x) =^ ^o-n^^^ 

n n 

in which a* is the complex conjugation of a„. Obviously f**{x) = f{x) and 
f{x)* = f*{x*). If a function satisfies /* = /, then it takes real values on the 
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real line. As mentioned above, e^^^ are shift operators e^^^fix) = f{x ^ and 
the Schrodinger equation 

'H(j)n{x) = 8{n)<P^{x), n = 0,l,2,..., (2.7) 

is a second order differential (oQM) or difference equation (dQM). 

Here are some explicit examples. From the ordinary QM. The prepotential w{x) 
parametrises the groundstate wavcfunction (f>a{x) which has no node and can be 
chosen real and positive <po{x) = e^^^K Then the above Hamiltonian takes the form 

H=p'' + U{x), U{x) = {dMx)f + d^w{x): 



(i) : Hermite 'w{x) = — |a;^, — oo < a; < oo, 

U{x) =x^ -\, T]{x) = X, (2.8) 

(ii) : Laguerre 'w{x) = —^x^ + glogx, 5 > 1, < a; < oo, 

U{x) = x^ + -l-2g, nix) = x^ (2.9) 

X'^ 

(iii) : Jacobi w(x) = glogsinx + /ilogcosx, 5 > 1, h > 1, < .t < 7r/2, 

I/(x) = ^^^^ + ^fc^-(5 + ft)2, ,,(:c)=cos2x. (2.10) 
sin'' a; cos^a; 

Prom the discrete QM with pure imaginary shifts (/3 = 7, < < 1): 

(iv): cont. Hahn V{x) = (ai + ix){a2 + ix), Iie{aj) > 0, —00 < a; < 00, 
7 = 1, V{x)=x, (2.11) 

t4 



ne^a^j > u, u < x 
7 = 1, {a^, Oj, Og, 04} = {ai, a2, as, 04} (as a set), r]{x) = x^ 



(v): Wilson V{x) = ^7] . -r-, R.e(aj) > 0, < x < 00, 



(2.12) 



n-=i(l-a,e«) 



(vi): Askey-Wilson F(x) = ^^-^5-^^^-— ^, K| < 1, < x < tt, 

7 = logg, {a^, ci3> CI4} = <i2, ct3, C14} (as a set), ry(x)=cosx. 

(2.13) 

From the discrete QM with real shifts (/3 = z): 

(vii): Hahn B{x) = (x + a)(iV - x), D(x) = x(6 + iV - x), a > 0, 6 > 0, 

n(x) = X, (2.14) 

(x + a)(x + 6)(x + c)(x + d) 



(viii): Racah -B(x) 
D(x) 



(2x + d)(2x + l + rf) 

(x + d — a) (x + d — 6) (x + d — c)x 
~ (2x - 1 + d)(2x + d) ' 

TV, a>6, 0<6<l + d, d>0, a>iV + d, 



J'^=i''a + 6 + c-d- 1, r7(x) = x(x + d), (2.15) 

/. ^ , r./ X (l-ag^)(l-6g^)(l-c9')(l-d9n 
(IX):, -Racah iJ(x) = (1 - dg4(l - ^ A ' ^''''^ 
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_ - (1 - a-Mg^)(l - 6-^dg^)(l - c~Mg"')(l - g"') 
(l-dg2a=-i)(i_rf^2x) 

c = q~^, a<b, 0<d<l, 0<a<q^d, qd<b<l, 

d<q-^, d = abcd-^q-^, 7?(a;) = (g"^ - 1)(1 - rfg^). (2.17) 

In the above list, the names like the Hermite, . . . , g-Racah are the names of the 
eigenpolynomials of the corresponding Hamiltonian in the form 

<f>n{x) = Mx)Pn{v{x)), (2.18) 

in which ri{x) is called the sinusoidal coordinate. It is also given in the same list. For 
the definitions and various properties of the hypergeometrie orthogonal polynomials 
and their g-versions in general, see (Koekoek & Swarttouw 1996) and (Odake & 
Sasaki 2008b, c) in connection with dQM. 

The groundstate wavefunction (l)o{x) is determined as a zero mode of A, 

AMx) = 0. (2.19) 

The similarity transformed Hamiltonian 7i in terms of the groundstate wavefunction 
4>o{x) has a much simpler form than the original Hamiltonian H: 

n =^ M^)~^ ° ^ ° M^) = £(v+{x){e^P - 1) + V-{x){e-^P - 1)) (2.20) 

^dw{x) d 
dx^ dx dx 

(v{x){e:'P -1) + V*{x){e-^'P -1)'^ pdQM • (2.21) 
(^B{x){l - e^) + D{x){l - e-^)^ pdQM 

It governs the differential (difference) equation of the polynomial eigenfunctions: 

nPnivix)) = £{n)Pn{n{x)). (2.22) 

Obviously, the square of the groundstate wavefunction (j)o{xY provides the orthog- 
onality weight function for the polynomials: 

oQM & pdQM I 'M^fPrMx))Pn{v{x))dx cx (5„„, (2.23) 

J Xi 

rdQM J2MxfPn{v{x))Pn{v{x))'xdnm- (2.24) 

Let me emphasise that the weight function, or cpQ (x) is determined as a solution of a 
first order differential (difference) equation, without recourse to a moment problem. 
This situation becomes crucially important when various deformations of orthog- 
onal polynomials are considered. The explicit forms of the squared groundstate 
wavefunction (weight function) ^o(a;)^ for the above examples in pure imaginary 
shifts dQM are: 

Mxf 



= < 
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r(ai + ix)T(a2 + ix)T{a'^ ~ ix)T{a2 — ix) : (iv) cont. Hahn 

4 

aj + ix)r{aj — ix) : (v) Wilson 

i=i 

4 

(e^" ; q)oo{e-'" ; q)ooY[{a,e''' ; q)^\a,e-'- ; : (vi) Askey- Wilson 

(2.25) 

For the real shifts dQM, the zero-mode equation A(j>o{x) — (j2.19p is a two term 
recurrence relation, which can be solved elementarily by using the boundary con- 
dition (l2Jl): 



: (vii) Hahn 



x\{N-xy. (6)jv 

(a,b,c,d)x 2x + d , .... ^ 

IVr-ri 1 \ 1^^ ^ : vhi Racah . (2.27) 

(l + d — a, 1-fa — 0, l-f-ct — c, l)a; a ^ ' 

1 - ^g'" . (ix^ o Racah 
, {a-^dq,b-^dq,c-^dq,q;q)j:d'' 1-a 

For the ordinary QM, the weight function is simply given by the prepotential 

(6) Intertwining Relations 

Let us denote by H^"' the original factorised Hamiltonian and by its partner 
Hamiltonian by changing the order of A'' and A: 

n^'^^^'A^A, n^'^^^'AAl (2.28) 

One simple and most important consequence of the factorisation is the intertwining 
relations: 

An^°^ = AA^A = hI^U, A^H^^^ = A^AA^ = n^°U\ (2.29) 

which is equally valid in the ordinary and the discrete QM. The pair of Hamiltonians 
'^['^1 and are essentially iso-spectral and their eigenfunctions {(j)n\x)} and 
{(/)L^'(a;)} are related by the Darboux-Crum transformations (Darboux 1882; Crum 
1955): 

n^"U^:!\x)=£{n)^^^\x), n = 0,l,..., (2.30) 

HW0W(a;) =f(n)0W(x), ti = 1,2,..., (2.31) 

0W(a;)=M°l(2:), 0l?l(2;) = ^0W(x), n=l,2,.... (2.32) 

The partner Hamiltonian T^I^l has the lowest eigenvalue £{1). If the groundstate 
energy is subtracted from the partner Hamiltonian 'H^^\ it is again positive 
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semi-definite and can be factorised in terms of now operators A^^^ and A^^^^: 

^[i]d|f_4[i]t^[i]+5(l). (2.33) 
By changing the orders of ^^^'^ and ^[^1, a new Hamiltonian "H^^l is defined: 

^[2] <^=^f^[iUWU£:(i). (2.34) 

These two Hamiltonians are intertwined by A^^^ and A^^^^: 

y^W (^[11 - £{!)) = ^[lUIilUIil = (7{[21 _ , (2.35) 

^[ilt(7^[2] _ f (1)) = ^Wt^m^Wt = (7^[i] _ £(i))^[ilt. (2.36) 

The iso-spectrahty of the two Hamiltonians 'HI^I and "H'^I and the relationship 
among their eigenfunctions follow as before: 

?^P]42](x)=£(n)42](x), n = 2,3,..., (2.37) 

</.Pl(x)=^W</.W(a:), 4^'(^) = g(^^''g(^) 4^'(^), n = 2,3,.... (2.38) 

This process can go on indefinitely by successively deleting the lowest lying energy 
level: 

n^''^<p'^Hx)^£{n)<f>^:\x), n = s,s + l,..., (2.39) 

2l[«-i]t 

4^1 (x) = ^[^-114-11 (x), 4-11 (x) = g(^)_g(,_i) 4^'(^). (2-40) 

for the ordinary QM (Crum 1955) and for the discrete QM as well (Gaillard & 
Matveev 2009; Odake & Sasaki 2009b). The determinant expressions of the eigen- 
functions (jM^ (x) are also known for the the ordinary QM (the Wronskian) and for 
the discrete QM as well (the Casoratian). 

By deleting a finite number of energy levels from the original Hamiltonian sys- 
tems and {(f)n\x)}, instead of the successive lowest lying levels, modification of 
Crum's theorem provides the essentially iso-spcctral modified Hamiltonian and 
its eigenfunctions {(f>n{x)}. The set of deleted energy levels V = {di, . . . ,di} must 
satisfy the conditions 

e 

'[[{m-dj)>0, m€Z+, (2.41) 

in order to guarantee the hermiticity (self-adjointness) of the modified Hamiltonian 
H and the non-singularity of the eigenfunctions 

n4)n{x) = £{n)4)n{x), neZ+\V. (2.42) 

Again the Wronskian (oQM) and Casoratian (pdQM) expressions of the eigenfunc- 
tions are known. For the ordinary QM, these results are known for some time (Krein 
1957; Adler 1994). For the dQM with pure imaginary shifts, the structure of the 
modified Crum's theory was clarified only very recently (Garcfa-Gutierrez, Odake 
& Sasaki 2010). Starting from an exactly solvable Hamiltonian, one can construct 
infinitely many variants of exactly solvable Hamiltonians and their eigenfunctions 
by Adler's and Garcia et al's methods. The resulting systems are, however, not 
shape invariant, even if the starting system is. For the dQM with real shifts, some 
partial results are reported in (Yermolayeva & Zhedanov 1999). 
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(c) Shape invariance 

Shape invariance (Gendenshtein 1983) is a sufficient condition for the exact solv- 
abihty in the Schrodinger picture. Combined with Crum's theorem (Crum 1955), 
or the factorisation method (Infold & Hull 1951) or the so-called supersymmetric 
quantum mechanics (Cooper, Khare & Sukhatme 1995), the totality of the discrete 
eigenvalues and the corresponding eigenfunctions can be easily obtained. It was 
shown by Odake and Sasaki that the concept of shape invariance works equally well 
in the discrete QM with the pure imaginary shifts (Odake & Sasaki 2005a, 2008c) 
as well as the real shifts (Odake & Sasaki 2008b), providing quantum mechanical 
explanation of the solvability of the Askey scheme of hypergeometric orthogonal 
polynomials in general. 

In many cases the Hamiltonian contains some parameter(s), A = (Ai, A2, . . .). 
Here we write parameter dependence explicitly, 'H(A), A{X), £{n; A), 0„(a;; A), etc, 
since it is the central issue. The shape invariance condition is 

A{X)A{Xy = nAiX + 6)^A{X + d) + f (1; A), (2.43) 

where k is a real positive parameter and S is the shift of the parameters. In other 
words H^^^ and H^^^ have the same shape, only the parameters are shifted by S. 
The energy spectrum and the excited state wavefunction are determined by the 
data of the groundstate wavefunction (poix; A) and the energy of the first excited 
state £(1; A) as follows: 

n-l 

£:(n;A) = ^k"£:(1;AW), aW=^A + s5, (2.44) 

s=0 

(t>n{x; A) oc yl(AM)ty^(AW)ty^(A[2l)t . . . A{X^'^-^^)Uq{x; A^). (2.45) 

The above formula for the eigenfunctions 0„(x;A) can be considered as the uni- 
versal Rodrigues type formula for the Askey scheme of hypergeometric polynomials 
and their g-analogues. For the explicit forms of the Rodrigues type formula for 
each polynomial, one only has to substitute the explicit forms of the operator A{X) 
and the groundstate wavefunction 4>o{x] A). For the nine explicit examples given in 
(|2.8p - (|2.17p . it is straightforward to verify the shape invariance conditions (|2.43p 
and the energy (|2.44p and the eigenfunction (|2.45l) formulas. For the ordinary QM, 
the above shape invariance condition (j2.43p becomes a relation among the prepo- 
tentials (k = 1): 

dwix;X)\'^ d^w(x;X) f dw(x; X + S)\'^ d^w(x]X + 6) m 

- ' ' ' ■ +t(l; A), (2.46) 



dx ) dx^ \ dx J dx'^ 

and the data for the three cases (i)-(iii) are 

(i): Hermite A = (null), £{n) = 2n, (2.47) 

(u): Laguerre A = c/, ^ = 1, £(n; A) = 4n, (2.48) 

(iii): Jacobi X={g,h), ^=(1,1), £{n; X) ^ 4n{n + g + h). (2.49) 

For the discrete QM with the pure imaginary shifts, the above shape invariance 
condition (|2.43p is rewritten as: 

V(x - q-, X)V*{x - q;X) = V{x; X + S)V*{x - ^7; A + d), (2.50) 
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V{x + i^-X) + V*{x- if ; A) = k{V{x- \ + 5) + V*{x\ X + 5))- £{!■ A). (2.51) 
The data for the three cases in pdQM (iv)-(vi) are (Odake & Sasaki 2005a, 2008c): 
(iv): cont. Hahn A=(ai,a2), ^ = (^,^), k = 1, 

£{n; A) = 4rt(ri + 6i — 1), bi '= ai + a2 + aj + Oj, 



(2.52) 

(v): Wilson A = (ai, a2, 03, 04) 5 = (i,i,i,i), k = 1, 



£■(«; A) = 4n(n + 61 — 1), 61 oi + 02 + 03 + 04, 



(2.53) 



(vi): Askey- Wilson = (ai, 02, 03, 04), 5 = (\,\,\,\), K = q^, 

8{n; A) = (q-" ~ 1)(1 - 649""^), 64 =^ 01020304. (2.54) 

For the discrete QM with the real shifts, the shape invariance (|2.43p is equivalent 
to the following set of two equations: 

B{x + 1; \)D{x + 1; A) B{x; A + S)D{x + 1; A + S), (2.55) 
B{x; A) + D{x + 1; A) = k{B{x; A + <5) + D{x; X + d))+ £{1; A). (2.56) 

The data for the three cases in rdQM (vii)-(ix) are (Odake & Sasaki 2008b): 

(vii) : Hahn A=(a,6,iV), ^ = (1,1,-1) k = 1, 

£{n;X) =4n{n + a + b~l), (2.57) 

(viii) : Racah X={a,b,d,N) ^ = (1,1,1,-1), k = 1, 

£{n;X)=4:n{n + d), (2.58) 

(ix) : (7- Racah q"^ — {a,b,d,q^^), ^ = (1,1,1,1), n — q^^ , 

£{n; A) = (g-" - 1)(1 - dV"^)- (2.59) 

It should be stressed that the size of the Hamiltonian {N in the finite case) decreases 
by one, since the lowest eigenstate is removed. 

(d) Solvability in the Heisenberg Picture 

All the Hamiltonian systems describing the hypergeometric orthogonal polyno- 
mials are exactly solvable in the Heisenberg picture, too (Odake & Sasaki 2006a, 
b). To be more precise, the Heisenberg operator of the sinusoidal operator ri{x) 

e'*"7/(a;)e^**^ (2.60) 

can be evaluated in a closed form. The sufficient condition for that is the closure 
relation 

[n, [H,v{x)] ] = r^ix) Roin) + [n, vi^)] Riiu) + R-iin). (2.61) 

Here Ri{y) are polynomials in y. It is easy to see that the cubic commutator 
\H, ["H, ["H, 77(0;)]]] = {a,<n-LYrj{x) is reduced to rj{x) and [H.,ri[x)\ with "H depending 
coefficients: 

{s.m)S{x) = [n,7^{x)]R^{n) + [H,[nMx)]]RiW 



Article submitted to Royal Society 



10 



= v{x) Rom)Ri{n) + [H, rtix)] {Ri{nf + Ro{n)) 

+ R-.i{H)Ri{H), (2.62) 

in which the definition {a.dTL)X =^ [TL,X] is used. It is trivial to see that all the 
higher commutators (ad'H)"77(x) can also be reduced to r]{x) and ['H,r]{x)] with TL 
depending coefficients. Thus we arrive at 



n=0 



in which the two "frequencies" a±{'H) are 



a±{n) = {Riin)±^/Ri{n)^ + ARo{n))/2, (2.64) 
a+{n) + a-{H) = Ri{H), a+{n)a-{H) = -Rq{H). (2.65) 

The annihilation and creation operators a^^^ are extracted from this exact Heisen- 

berg operator solution: 



e 



r]{x)e-'''^ = „(+)g»a+(w)t ^ ^(-)g»a_(w)t _ (2.66) 

''=^'±([^,77(0;)] - (ry(x) + i?_i(^)i?o(^)-')a^(^)) (a+(^) - Q-(^))"' 

= ±(a+(?^) - a-{H))~'{[n,v{x)]+a±{n){v{x)+R-i{H)Ro{H)-^)). 

(2.67) 

The energy spectrum is determined by the over-determined recursion relations £{n+ 
1) = S[n) + a+{£{n)) and £{n - 1) = £{n) + a-{£{n)) with £{0) = 0, and 
the excited state wavefunctions {(^„(x)} are obtained by successive action of the 
creation operator a(+^ on the groimdstate wavefunction (l>o{x). This is the exact 
solvability in the Heisenberg picture. 

The data for the three cases in oQM (i)-(iii) are 

(i) :Hermite Ro{y) = 4, i?i(y) = i?_i(y) = 0, (2.68) 

(ii) : Laguerre Ro{y) = 16, Ri{y)=0, R-i{y) = -8{y + 2g + 1), (2.69) 

(iii) : Jacobi Ro{y) = 16(y + {g + hf), Ri{y) = 0, 

R-i{y) = 16{g-h){g + h-l). (2.70) 

The data for the three cases in pdQM (iv)-(vi) are (Odake & Sasaki 2005a, 2008c): 

(iv) : cont. Hahn Ro{y) = 4y + 4Re(ai + a2)(Re(ai + 02) - l), Ri{y) = 2, 

R-i{y) = 2Im(ai + a2)y + 4(Re(ai + 02) - l)lm(aia2), (2.71) 

(v) : Wilson R^{y) = 4y + 61(61 - 2), Ri{y) = 2, 

R-i{y) = -22/2 + (61 - 262)2/ + (2 - 61)63, 
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dcf ^ dcf ^ ajQkai, (2.72) 

l<i<'=<4 l<j<fc</<4 

(vi): Askey- Wilson Ri{y) = {q-^ - q^ fy', y' =' y + 1 + q^^bi, 

R.i{y) = ^^{q-i ~ q^)'{{h + q-%)y' - (1 + q-'Kb^ + q-^b.h^)), 

4 

6i=^aj, 63= ^ ajakai. (2.73) 

i=i i<j<fc<i<4 

The data for the three cases in rdQM (vii)-(ix) are (Odake & Sasaki 2008b): 

(vii) : Hahn i?o(y) = 4y + (a + 6 - 2)(a + 6), i?i(y) = 2, 

i?_i(y) = -(2Ar - a + 6)2; - a(a + 6 - 2)iV, (2.74) 

(viii) : Racah i?o(y) = 4y + - 1, Ri{y) = 2, 

i?_i(y) = 2y2 + (2(a6 + 6c + ca) - (1 + d)(l + d))y + abc{d - 1), (2.75) 

(ix) : <7-Racah Roiy) ^ {q-^ ~ q^^ {y'^ - (q-^- + q^ fd) , y' = y + 1 + d, 

Ri{y) = (g-^ -g^)'y', 

R-iiy) = (g^^ - g^)2|^(i + cf)?;'2 - (a + 6 + c + d + J+ (a5 + 6c + ca)g"^)?/' 
+ ((l-a)(l-6)(l-c)(l-dVi) 

+ (a + 5 + c - 1 - + (a6 + 6c + ca)q-^)) (1 + J)) . (2.76) 

(e) Closure Relation 

The dwa/ closure relation has the same forms as the closure relation (I2.6ip with 
the roles of the Hamiltonian H and the sinusoidal coordinate r]{x) exchanged: 

[77, [77, n]] = n R^r\v) + bh n] r^,^^\v) + R^-T'iv). (2.77) 

in which 

(77(2; - il3) - r^ix)) + {ii{x + ip) ~ 77(x)) , (2.78) 

^{r^{x - i(3) - 77(2;)) (77(2; + i(3) - r^ix)) , (2.79) 

e{V+ix) + V^ix))R^,^^\rj{x)). (2.80) 

The dual closure relation is the characteristic feature shared by all the 'Hamilto- 
nians' Ti which map a polynomial in 77(2:) into another. Therefore its dynamical 
contents are not so constraining as the closure relation, except for the real shifts 
(the discrete variable) exactly solvable {L — 2) case, where the closure relation 
and the dual closure relations are on the same footing and they form a dynamical 
symmetry algebra which is sometimes called the Askey- Wilson algebra (Zhedanov 
1992; Granovskii et al 1992; Terwilliger 2004; Koornwinder 2008; Odake & Sasaki 
2008b). 



R''_r\vi^)) = 
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3. Unified Theory of Exact and Quasi-Exact Solvability 

In the previous section various examples of exactly solvable systems arc explored 
and shown that the two sufficient conditions for exact solvability, the shape invari- 
ance and the closure relation, are satisfied by them. However, these two sufficient 
conditions do not tell how to build exactly solvable models. In this section I present 
a simple theory of constructing exactly solvable Hamiltonians in discrete QM. It 
covers all the known examples of exactly solvable discrete QM with both pure imag- 
inary and real shifts and it predicts several new ones to be explored. The theory is 
general enough to generate quasi- exactly solvable Hamiltonians in the same manner. 
The quasi-exact solvability means, in contrast to the exact solvability, that only a 
finite number of energy eigenvalues and the corresponding eigenfunctions can be ob- 
tained exactly. Many examples are known in the ordinary QM (Morozov et al 1990; 
Ushveridze 1994) but only a few are known in dQM (Wiegmann & Zabrodin 1994, 
1995) . This unified theory also incorporates the known examples of quasi-exactly 
solvable Hamiltonians (Odake & Sasaki 2007; Sasaki 2007, 2008). A new type of 
quasi-exactly solvable Hamiltonians is constructed. The present approach reveals 
the common structure underlying the exactly and quasi-exactly solvable theories. 
This section is a brief introduction to a recent work by Odake and Sasaki 2009a. 

In the following I will take the similarity transformed Hamiltonian Ti. (|2.20p 
instead of "H as the starting point. That is, I reverse the argument and construct 
directly the 'Hamiltonian' 71 (j2.2ip based on the sinusoidal coordinate r](x). The 
general strategy is to construct the 'Hamiltonian' T-l in such a way that it maps a 
polynomial in ?/(x) into another: 

C Vn+L-2 C Voo ■ (3.1) 

Here V„ (n e Z>o) is defined by 

Vn Span[l, 77(x), . . . , vixT] , Voo Hm V„. (3.2) 



(a) potential functions 

The general form of the 'Hamiltonian' H mapping a polynomial in ri{x) into 
another is achieved by the following form of the potential functions V± (x) : 

(-^^ ^ y±{x) 3-) 

{rj{x =F ip) ~ ri{x)) (ri{x ^ i/S) — rj{x ± ifi)) ' 
V±{x)= Vk,i^{x)S{xTiP)\ (3.4) 



k + l<L 



where L is a natural number indicating the degree of ri[x) in V±{x) and vj^^i are 
real constants, with the constraint ^ ; 7^ 0. It is important that the same Vk,i 

k+l=L 

appears in both V±{x). The 'Hamiltonian' T-L with the above V±{x) maps a degree 
n polynomial in ri{x) to a degree n -\- L — 2 polynomial. 

The essential part of the formula (13. 3p is the denominators. They have the same 
form as the generic formula for the coefficients of the three term recurrence relations 
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of the orthogonal polynomials, (4.52) and (4.53) in (Odake & Sasaki 2008b). The 
translation rules are the duality correspondence itself, (3.14)-(3.18) in (Odake & 
Sasaki 2008b): 

£{n) ri{x), ~An^V+{x), -C„^y_(a;), 
a+[£{n)) ^ rj{x — i/3) — ri{x), (£(?i)) — > 77(x + i/3) — 77(2;). (3.5) 

Some of the parameters Vk,i in (13.41) are redundant. It is sufficient to keep Vk^i 
with ^ = 0,1. The remaining 2L + 1 parameters Vk,i {k + I < L, I = 0,1) are 
independent, with one of which corresponds to the overall normalisation of the 
Hamiltonian. 

The^L = 2 case is exactly solvable. Since the Hamiltonian of the polynomial 
space "H is an upper triangular matrix, its eigenvalues and eigenvectors are eas- 
ily obtained explicitly. For the solutions of a full quantum mechanical problem, 
however, one needs the square-integrable groundstate wavefunction (t>o{x) (j2.19p . 
which is essential for the existence of the Hamiltonian H and the verification of 
its hermiticity. These conditions would usually restrict the ranges of the param- 
eters wo,Oj ■ • • J ^2.0- It is easy to verify that the explicit examples of the potential 
functions V{x), V*{x) and B{x), D{x) in (j2.1ip - (j2.17p are simply reproduced by 
proper choices of the parameters {wfe,;}. It should be stressed that the above form 
of the potential function (I3.3p - (|3.4I) provides a unified proof of the shape invariance 
relation (|2.43p . the closure relation (I2.6ip and the dual closure relation (|2.77p in 
the % scheme. 

The higher L > 3 cases are obviously non-solvable. Among them, the tame 
non-solvability of L = 3 and 4 can be made quasi-exactly solvable (QES) by adding 
suitable compensation terms. This is a simple generalisation of the method of Sasaki 
& Takasaki 2001 for multi-particle QES in the ordinary QM. For a given positive 
integer M, let us try to find a QES 'Hamiltonian' or more precisely its modifi- 
cation "H', having an invariant polynomial subspace Vm'- 

U'Vm C Vm- (3.6) 

For L ^ 3, TL' is defined by adding one single compensation term of degree one 

n'^^'n-eoiMUx), (3.7) 

and we have achieved the quasi-exact solvability TL'Vm Q Vm- Known discrete QES 
examples belong to this class (Sasaki 2007, 2008). 

For L = 4 case, H' is defined by adding a linear and a quadratic in ri{x) com- 
pensation terms to the Hamiltonian H: 

n' - eoiM)vixf - ei{M)r){x), (3.8) 
and n'r){xy^ G Vm- This type of QES theory IS new. 
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